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Abstract: We study the existence of non-collision periodic solutions with Newtonian 
| potentials for the following planar restricted 4-body problems: Assume that the given 

CN ■ positive masses 7711,7712,7713 in a Lagrange configuration move in circular obits around 

their center of masses, the sufficiently small mass moves around some body. Using 
variational minimizing methods, we prove the existence of minimizers for the Lagrangian 
action on anti-T/2 symmetric loop spaces. Moreover, we prove the minimizers are non- 
collision periodic solutions with some fixed wingding numbers. 
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1 Introduction and Main Results 



In this paper, we study the planar circular restricted 4-body problems with Newtonian 
■ potentials. Suppose points of positive masses 7711,7712,772-3 move in a plane of their circular orbits 

Qi{t), ?2(i)j Qs{t) and the center of masses is at the origin; suppose the sufficiently small mass 
I point does not influence the motion of 777 1 , 777,2 , 7773 , and moves in the plane for the given masses 

7771, 777 2 , 777 3 . 

. It is well-known that qi(t),q2(t),q3(t) satisfy the Newtonian equations: 

CO 

dU 

miqi = —, 7 = 1,2,3, (1.1) 
% 

^ ■ where 

u= y, , m : m V (i-2) 

Without loss of generality, we assume that there exists 61,82,63 £ [0, 27r) such that the planar 
circular orbits are 

qi (t) =rie >FT^t e V=T(? lj q2{t)=r2e ^i^t e V^ie^ q ^ = rg< y=T^t e V=Tfl Sj (L3) 

where the radius r\,r2,r 3 are positive constants depending on ?77j(7 = 1,2,3) and T (see Lemma 
2.6). We also assume that 

m iqi {t) + m 2 q2(t) + m 3 q 3 (t) = (1.4) 

and 

\qi-qj\ = i, i<£^j<3, (1.5) 

where the constant I > depends on mj(i = 1, 2, 3) and T (see Lemma 2.5). 
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The orbit q(t) £ R 2 for sufficiently small mass is governed by the gravitational forces of 
mi,m.2,m3 and therefore it satisfies the following equation 



3 

£ 

i=i 



miqi - q) 



(1.6) 



For iV-body problems, there are many papers concerned with the periodic solutions by using 
variational methods, see [1-9,13-16,18] and the references therein. In [3], Chenciner- Montgomery 
proved the existence of the remarkable figure- "8" type periodic solution for planar Newtonian 3- 
body problems with equal masses. Marchal[6] studied the fixed end problem for Newtonian n-body 
problems and proved the minimizer for the Lagrangian action has no interior collision. Especially, 
in [8], Simo used computer to discover many new periodic solutions for Newtonian n-body problems. 
Zhang-Zhou[13-15] decomposed the Lagrangian action for n-body problems into some sum for two 
body problems and [14,15] avoid collisions by comparing the lower bound for the Lagrangian action 
on the symmetry collision orbits and the upper bound for the Lagrangian action on test orbits in 
some cases. 

Motivated by the above works, we use variational methods to study the circular restricted 
3+1-body problem with some fixed wingding numbers and some masses. 

For the readers' conveniences, we recall the definition of the winding number, which can be 
found in many books on the classical differential geometry. 

Definition 1.1 Let T : x(t),t 6 [a, b] be an given oriented continuous closed curve, and p be 
a point of the plane not on the curve. Then, the mapping tp : T — > S 1 , given by 

y(g(t))= , X fl~ P , , t€[a,b] 
\x{t) — p\ 

is defined to be the position mapping of the curve T relative to p, when the point on T goes around 
the curve once, its image point <p{x(t)) will go around S 1 a number of times, this number is called 
the winding number of the curve T relative to p, and we denote it by deg(F,p). If p is the origin, 
we write degT. 
Define 

W 1 ' 2 {R/TZ,R 2 ) = lx(t) x(t),x{t) e L 2 (R,R 2 ), x(t + T) =x(t)\. 
The norm of W 1 ' 2 (R/TZ, R 2 ) is 



r f T 


i 


- fT 1 


\ \x\ 2 dt 


5 + 


/ \x\ 2 dt 


L Jo 




-Jo J 



(1.7) 



The functional corresponding to the equation (|1.6|) is 



where 



A. 



q G W l > 2 {R/TZ,R z ) 



T 



and 



A + = K G W h2 (R/TZ,R 2 ] 
Our main results are the following: 



i=l 



mi 



dt, q G A±, 



lit + I) = deg(q - qi) = -1, 

q(t)^ qi (t), Vie [0,T],i = 1,2,3 



q{t + |) = ~q{t), deg(q - q{) = 1, 
q(t)^ qi (t), ViG [0,T],i = 1,2,3 
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Theorem 1.1 Let T = 1, for the values of 7711,7712,7713 given in Table 1 with M = 1, the 
minimizer of f(q) on the closure A_ of A_ is a non-collision 1-periodic solution of fjl .6[) : for the 
values of mi, 7712, ms given in Table 2 with 777,1 = 7712 = 7773 = 1, the minimizer of f(q) on A_ is a 
non-collision 1-periodic solution of (|l,6p . 



Remark 1 In proving Theorem 1, we need to use test functions. We find that if the test 
functions are circular orbits, we can not get the desired results on A_. Therefore, we select elliptic 
orbits as test functions. 

Theorem 1.2 Let T = 1, for the values of 7771,7712,7713 given in Table 3 with M = 1, the 
minimizer of f(q) on the closure A + of A + is a non-collision 1-periodic solution of fjl .6[) : for the 
values of 777,1,7772,777,3 given in Table 4 with 777-1 = 1^2 = ^3 = 1, the minimizer of f(q) on A + is a 
non-collision 1-periodic solution of 



Remark 2 When we take elliptic orbits as test functions, we find that the biggest symmetric 
space is the anti-T/2 symmetric loop space if the wingding number 77 is odd(n = ±1, ±3, • • • ); we 
can not find suitable symmetric space if the wingding number is even. When the wingding number 
77 ^ ±1 and we take circular orbits as test functions, we find that the biggest symmetric space is 

A= f ll2(i?/ 2) q(t + ] ^)=R(\n-l\)q(t), deg( q - qi )=n, \ 

V q(t)^qi(t), Vie [0,T],7 = 1,2,3 f 



where 



cosj^rr sin, 271 



^ |n " 1|)= ( S »H C0S H l )^SO(2) 



is a counter-clockwise rotation of angle i^rn i n R 2 - But the Lagrangian actions on the circular test 
orbits are bigger than the lower bound for the Lagrangian actions on collision symmetric orbits. 
Hence we consider the anti-T/2 symmetric loop spaces A±. 

2 Preliminaries 

In this section, we will list some basic Lemmas and inequality for proving our Theorems 1.1 and 
1.2. 

Lemma 2 . 1 (Tonelli [Tj . [TT] ) Let X be a reflexive Banach space, S be a weakly closed subset 
of X, f : S — > RL) {+00}. If / ^ +00 is weakly lower semi-continuous and coercive(/(x) — > +00 as 
||x|| — > +00), then / attains its infimum on S. 

Lemma 2.2(Poincare-Wirtinger Inequality [TO]) Let q G W 1 ' 2 (R/TZ, R K ) and q(t)dt = 0, 

then 



dt. 



Lemma 2.3(Palais's Symmetry Principle([12])) Let a be an orthogonal representation of a 
finite or compact group G, H be a real Hilbert space, / : H — >■ R satisfies f(a ■ x) = f(x),\/a € 
G,Vx e H. 

Set F = {x £ H\a ■ x = x, Vc £ G}. Then the critical point of / in F is also a critical point 
of / in H. 

Remark 2.1 By Palais's Symmetry Principle and the perturbation invariance for wingding 
numbers, we know that the critical point of f{q) in A± is a periodic solution of Newtonian equation 

<HH). 
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Lemma 2.4 

(1) (Gordon's Theorem [17]) Let x £ W 1 ' 2 ([t 1 ,t 2 ], R K ) and x(h) = x(t 2 ) = 0. Then for any 
a > 0, we have 

f\\\*\ 2 + A)* > ^(2vr) 2 / 3 a 2 / 3 (t 2 - t^Vs. 

Jti 2 M 2 

(2) (Long-Zhang[T8]) Let x € W^ 2 (R/TZ, R K ), xdt = 0, then for any a > 0, we have 

/" T (i|i| a + ^-)dt > -(2^) 2 / 3 a 2 / 3 T 1 /3. 

jo 2 fI 2 



Lemma 2.5 Let M = mi + m 2 + m3, we have Z = y , 
Proof. It follows from JTI]) and JT2J) that 



qi = m 2 -, pr + ms- R . (2.1) 

92-91 93-91 



Then by ([T5 ]) -(fI3 ]) . we obtain 



which implies 



47T 2 1 . 

--y2"9l = p("T-292 + "1393 - ™29l ~ ™39l, 

= ^(-miQi - m 2 gi - m 3 gi), 



,3 _ MT 2 



4^ 2 ' 



(2.2) 



(2.3) 



that is, 



3 MT 2 

Lemma 2.6 The radius r\, r 2 , r% of the planar circular orbits for the masses mi, m 2 , m 3 are 



Vm 2 +m 2 m 3 +m 2 
ri = M 



v/m? + mim 3 + ml , 
r2 = M 

-v/m? + mim 2 + m| , 

rs = — = -I. 

A M 

Proof. Choose the geometrical center of the initial configuration [q\ (0) , q 2 (0) , g 3 (0) ) as the ori- 
gin of the coordinate (x,y). Without loss of generality, by (|1.5p . we suppose the location coordinates 
of 9i(0),g 2 (0),g 3 (0) are A 1 (^-,0),A 2 (-^,^),A 3 (-^,-^). Then we can get the coordinate of 

the center of masses mi, m 2 , m 3 is C(— ^ , 2 M 2 ). To make sure the Assumption 

(|1.4p holds, we introduce the new coordinate 



Y _ _ &m 1 l-&m 2 l-^m 3 l 
A — X M , 

7712 7 m ^ 



y = y-^ — 2- 



IF 
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Hence in the new coordinate (X,Y), the location coordinates of qi(0), 92(0), 93(0) are A\{ 



A 2 (- 



&mil ^l+m 3 U . , &mil ^l+m 2 l 



M 



,4 a 



M /' M ' M J'^\ M > M 

at the origin 0(0,0). Then compared with (|1.3p . we have 



) and the center of masses mi, 7712,7713 is 



and 



sin 0i 



sin 02 



sin 03 



n = |AiO| 

r 2 = |A 2 0| 
r 3 = |A 3 0| 

771 2 + "13 



■y/ 7772 + m 2 m 3 + m 3 

M 

1/ mf + minis + 7713 
M 



\/ m^ + 77117712 + 



1y/m\ + 71727713 + 77I3 
Till + 27713 

2y/m\ + 77117713 + 7173 

771i + 27112 

"2V^+mim 2 +m^ 



M 



COS 01 



COS 02 



COS 03 



\/3(m 2 + 



7713 ) 



2^/m\ + 717-27713 + 77I3 

\/3mi 
2^Jm\ + mim?, + mg 



\/3mi 



2 Vm? + mim 2 + m^ 



□ 



(2.5) 

(2.6) 
(2.7) 

(2.8) 
(2.9) 
(2.10) 



3 Proof of Theorems 



In order to get Theorems, we need two steps to complete the proof. 

Step 1: We will establish the existence of variational minimizers of f(q) in (|1.8p on A±. 

Lemma 3.1 f(q) in (|1.8p attains its infimum on A±. 

Proof. By using Lemma 2.2, for \/q G A±, we can get that the equivalent norm of (|1.7|) in A± 



is 



\q\ 2 dt 



(3.1) 



Hence by the definition of f(q), f is coercive on A±. Next, we claim that / is weakly lower semi- 
continuous on A±. In fact, for \/q k G A±, if q k — q weakly, by compact embedding theorem, we 
have the uniformly convergence: 



max \q (t) - q(t)\ -> 0, k — > oo, 
0<t<T 



which implies 



T 3 
i=l |y 



■mi 



-dt 



T 3 

E 

i=l 



777; 



(3.2) 



(3.3) 



It is well-known that the norm and its square are weakly lower semi-continuous. Therefore, com- 
bined with (|3.3p . we obtain 

lirniaf /(?*)>/(?), (3.4) 

k— >oc 

that is, / is weakly lower semi-continuous on A±. By Lemma 2.1, we can get that f{q) in (|1.8p 
attains its infimum on A±. □ 

Step 2: We will prove the variational minimizers in Lemma 3.1 is the noncollision T-period 
solution of (|1.6p . 

For any collision generalized solution q, we can estimate the lower bound for the value of 
Lagrangian action functional. 
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Lemma 3.2 For dA± = {q G W 1 > 2 (R/TZ, R 2 )\q{t + f) 
[0,T] s.t. g^i^i) = g(*j±)}> we have 



-?(*), 31 < »± < 3,t.± e 



inf f(q) > -(2^) 2 / 3 CM^ 1 /3 T i/3 A ^ 
g69A± 2 



where 



2 f mi + m2 + m3 _ ^( mim2 + mi m 3 + m 2 m 3 ), 
C = min < 23?7i2 + mi + 7713 — g^y(?7ii?7i2 + mim.3 + i^mz), 
k 2l m3 +m 1 +m 2 - ! l 7 (m ira2 + m 1 m3+ m2m 3) 



Proof. It follows from (|1.4p that 



which implies 



mift = 0, 



i=l 



(3.5) 



i=i 



^rn i (|(?| 2 + |g l | 2 -2(g,g i ) 
t=l 

3 

M\q\ 2 + ^m i \q i \ 2 -2^q,^2m iqq 
i=l 
3 

M\q\ 2 + ^m\qi\ 2 . 



i=l 



(3.6) 



Therefore 



Hence 



/(<?) 



1 3 







1 



• 1 2 



i=l 



171; 



(It 



^ m \^-^ + \q- qi 



1 f T 3 rl 



dt 



2M 



T 3 



1=1 



(3.7) 



(3.8) 



If g 6 A_ is a collision generalized solution, then there exists t-- € [0, T] and 1 < i < 3 such that 
= <?;-(%)• Since + f ) = -9,(0, we obtain q(t^ +f ) = + f ), V0 < < 2. 

So, by (1) of Lemma 2.4, we get 



M 



m- 



l o 12 



■ ,2 M 



dt 



-m- 



M J o 7 L2 



1,. . |2 M 



;l9-? in -l ' + 



19 - ?i 



> £(27r) 2 / 3 2 2 / 3 m 4 -Af- 1 /3 T V3. 



(3.9) 



For noncollision pair q,qi(i 7^ z ), we have J Q q(t)dt = 0, J Q qi(t)dt = 0. Therefore J Q T (<?(£) 
qi(t))dt = 0. Hence by (2) of Lemma 2.4, we can get 



Til I>[^-*l 2 + 



M 



k - % 



> |(2^) 2 / 3 (A/ - JM-^T 1 / 3 . (3.10) 
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For the other term of /, using the expression for the orbits 9i,9 2 ,9 3 as in (|1.3p . Lemma 2.5 and 
Lemma 2.6, we obtain 

I fT 3 i 
~2M Y, mi ^\ 2dt = " 2 ( 2 ^) 2/3 (™i m 2 + mim 3 + m 2 m 3 )M~ 4 / 3 T 1 / 3 . (3.11) 



i=l 



Therefore, it follows from (|3.9p - (|3.1ip that 



inf f(q) > -(2^) 2 / 3 CM~ 1 / 3 r 1 /3 4 ^ 
qGdA_ 2 



(3.12) 



where 



f 2 i 

2smi + m 2 + m 3 - 3^(777,1777^ + mim 3 + m 2 m 3 ), 

C = min ^ 23rri2 + mi + 777.3 — 357(^1^2 + minis + 777, 2 77i 3 ), 
2 1 
2 3 777,3 + mi + m 2 - 3^(771,17712 + miffl 3 + m 2 m 3 ) 

Similarly, if q € A + is a collision generalized solution, we have 



inf /(g) > -(2vr) 2 / 3 CM- 1 / 3 T 1 / 3 4^, □ 

g6<9A+ 2 



(3.13) 



Proof of Theorem 1.1 In order to get Theorem 1.1, we are going to find a test loop q € A_ 
such that f(q) < d 2 . Then the minimizer of / on A_ must be a noncollision solution if ci 2 < d\. 

Let a > 0, 6 > 0, 9 £ [0, 2tt) and 



q — qi = I 
Hence 



2vr \ , . / 2tt 
— i + 0j,&sm(^- — t + 



(3.14) 



9-92 = 9-91+91-92 

= (9i - 92) + (9- 9i) 



ri cos ( — t + 9i) — r 2 cos f — t + 2 ) + a cos ( — — t + 9 ) , r\ sin ( — i + 9\ 



\ T 



\ T 



2vr 



/2vr 



-r 2 sin + 6> 2 ) + 6sin ( - -^t + 



(3.15) 



9-93 



( Tpi + #1) - ^3 cos + 6> 3 ) + a cos ( - + 6>) , n sin + 6>i 



2yT 

-r 3 sin ( — i + 6> 3 ) + 6 sin 



2vr 
~T 



t + < 



(3.16) 



It is easy to see that q G A_ and 



19 - 91 



2vr\2 



a 2 + b 2 a 2 -b 2 Mir 

cos(—t - 29 

2 2 V T 



(3.17) 



19 - 9i I 



a 2 + b 2 



a 2 -b 2 /47T 



(3.18) 
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\Q-Q2\ 



2tt\2 f g 2 + fr 2 
TV 1 2 



-y-^cos^t - 20^ + r 2 + rf - 2rir 2 cos(0 2 



r /47T \ /47T 

-(a + 6) ^ncos^— t + 0\ - 9) - r 2 cos(^— t + 2 
+(a - b) [ncos(0i + 0) - r 2 cos(9 2 + 0)] 



\Q ~ 12 



a 2 + b 2 a 2 - b 2 



+ 



4tt 



-cos 



2 

47T 



(^t -29)+ri + r z 2 - 2 ri r 2 cos{9 2 - 9 X ) 



I -Tpt + 9 1 -9j- r 2 cos\^—t + 9 2 



+ (a + b) 

+(a - 6) [ncos((9i + 9)- r 2 cos{9 2 + 0)] 



47T 



.. . |2 /2vr\2 fa 2 + 6 2 

l*-<»l = [jr) |— 2- 



47T 



-(a + 6) 

+(a - 6) [r l cos(9 1 + 9)- r 3 cos(9 3 + 0)] 1 



47T 

r 3 cos(— t + 3 -0 



9 - 93 



a 2 + b 2 | a 2 



4"7T 



-cos 



t - 29 + r 2 + r 2 - 2rir 3 cos(0 2 - 0i) 



+ (a + 6) 



4-7T 

ricosl — t + 6'i-6' 



47T 

r 3 cosl — « + 03 



T / J \ T 

+(a - 6) [ricos(0i + 0) - r 3 cos(0 3 + 0)] 



m 2 - (|) 2 r ;, = (|) 2 -i w = (f)V§. 



Therefore by (1X171) - (fc^) . we get 

(■T 3 

M 



/(«) 



M 



\q-Qi 



dt 



2M 



T 3 



t=i 



mi | dil dt 



2ir ( a +b m 2 + 777,3 — m i 2 2m 2 r 2 cos(9 2 — 9\) + 2?n3r3cos(03 — 0i" 
j 1 — r^ 



+ 



T I 2 

m 2 (a - 6) 



M 



M 



M 



[r 1 cos{9 1 + 0) - r 2 cos(0 2 + 0)] + ^ [r lC os(0i + 0) 



-r 3 cos(0 3 + 0)]| + 777-1 



T 



M 

a 2 + 6 2 a 2 -6 2 /4vr n/1 

— 1 cos I — 1-29 

2 2 \T 



dt 



A f T (a 2 + b 2 a 2 -b 2 



i=2 



2 
47T 



47T 



+ (a + b) [r lC os^t + 



- — cos \^—t - 29 ) + r 1 + r\ - 2r x riCos{vi 



'4?T 

riCos[ —t + 9i- 9 



+(a - b) [r 1 cos(9 1 + 0) - r iC os(0i + 0)] 
d 2 (a,b,9). 



dt 



In order to estimate d,2, we have computed the numerical values of di = f(q) over some selected 
test loops. The computation of the integral that appears in (|3.24p has been done using the function 
{quad} of Mathematica 7.1 with an error less than 10 -6 . Let T = 1, the results of the numerical 
explorations are given in Table 1 with M = 1 and Table 2 with mi = mi = m% = 1. 

Table 1: Parameters for test loops for Theoerm 1.1 



a 


b 


e 


ml 


m2 


Hl3 


dl 


d2 


0.13 


0.49 


7T 

20 


0.29 


0.42 


0.29 


5.419669 


5.417862 


0.15 


0.49 


7T 

20 


0.29 


0.41 


0.30 


5.417626 


5.416591 


0.15 


0.49 


7T 

20 


0.29 


0.42 


0.29 


5.419669 


5.413794 


0.15 


0.49 


7T 

20 


0.30 


0.35 


0.35 


5.441499 


5.436767 


0.15 


0.51 


7T 

20 


0.30 


0.36 


0.34 


5.441669 


5.437985 


0.15 


0.51 


7T 

20 


0.30 


0.37 


0.33 


5.442180 


5.433615 


0.15 


0.51 


7T 

20 


0.30 


0.38 


0.32 


5.443031 


5.429587 


0.15 


0.51 


7T 

20 


0.30 


0.39 


0.31 


5.444223 


5.425898 


0.15 


0.51 


7T 

20 


0.30 


0.40 


0.30 


5.445755 


5.422550 


0.15 


0.53 


7T 

20 


0.31 


0.35 


0.34 


5.470820 


5.467576 


0.15 


0.53 


7T 

20 


0.31 


0.36 


0.33 


5.471160 


5.462971 


0.15 


0.53 


7T 

20 


0.31 


0.37 


0.32 


5.471841 


5.458707 


0.15 


0.53 


7T 

20 


0.31 


0.38 


0.31 


5.472863 


5.454784 


0.17 


0.45 


7T 

20 


0.32 


0.32 


0.36 


5.500992 


5.488608 


0.17 


0.47 


7T 

20 


0.32 


0.33 


0.35 


5.500481 


5.454518 


0.17 


0.47 


7T 

20 


0.32 


0.34 


0.34 


5.500311 


5.449987 


0.17 


0.47 


7T 

20 


0.33 


0.34 


0.33 


5.530142 


5.444254 


0.45 


0.15 


7T 


0.33 


0.31 


0.36 


5.471160 


5.456006 


0.45 


0.15 


7T 


0.33 


0.32 


0.35 


5.500481 


5.455325 


0.45 


0.15 


7T 


0.33 


0.33 


0.34 


5.530142 


5.454984 


0.47 


0.13 


7T 


0.34 


0.30 


0.36 


5.441669 


5.439671 


0.47 


0.13 


7T 


0.34 


0.31 


0.35 


5.470820 


5.438820 


0.47 


0.13 


7T 


0.34 


0.32 


0.34 


5.500311 


5.438309 


0.47 


0.15 


7T 


0.35 


0.30 


0.35 


5.441499 


5.417900 


0.49 


0.15 


7T 


0.36 


0.29 


0.35 


5.412519 


5.411552 


0.49 


0.15 


7T 


0.36 


0.32 


0.32 


5.500992 


5.410020 


0.49 


0.15 


7T 


0.37 


0.29 


0.34 


5.412859 


5.411962 


0.49 


0.15 


7T 


0.37 


0.30 


0.33 


5.442180 


5.411281 


0.49 


0.15 


7T 


0.37 


0.31 


0.32 


5.471841 


5.410940 


0.49 


0.15 


7T 


0.38 


0.29 


0.33 


5.413540 


5.412712 


0.49 


0.15 


7T 


0.38 


0.30 


0.32 


5.443031 


5.412201 


0.49 


0.15 


7T 


0.38 


0.31 


0.31 


5.472863 


5.412031 


0.49 


0.15 


7T 


0.39 


0.29 


0.32 


5.414562 


5.413803 


0.49 


0.15 


7T 


0.39 


0.30 


0.31 


5.444223 


5.413462 


0.49 


0.17 


7T 


0.40 


0.29 


0.31 


5.415924 


5.415807 


0.49 


0.17 


7T 


0.40 


0.30 


0.30 


5.445755 


5.415637 


0.49 


0.17 


7T 


0.41 


0.30 


0.29 


5.417626 


5.416078 


0.49 


0.17 


7T 


0.42 


0.29 


0.29 


5.419669 


5.416689 
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Table 2: Parameters for test loops for Theoerm 1.1 



a 


b 


e 


ml 


m2 


m3 


dl 


d2 


0.15 


0.67 


7T 

,-su 


1.00 


1.00 


1.00 


11.523843 


11.505860 


0.15 


0.67 


7T 
,11 J 


1.00 


1.00 


1.00 


11.523843 


11.505860 


0.15 


0.69 


7T 

,-su 


1.00 


1.00 


1.00 


11.523843 


11.444212 


0.17 


0.65 


7T 
rfO 


1.00 


1.00 


1.00 


11.523843 


11.493238 


0.17 


0.67 


7T 

20 


1.00 


1.00 


1.00 


11.523843 


11.452135 


0.17 


0.69 


7T 

20 


1.00 


1.00 


1.00 


11.523843 


11.400124 


0.19 


0.63 


7T 
rfO 


1.00 


1.00 


1.00 


11.523843 


11.519350 


0.19 


0.65 


7T 

on 
20 


1.00 


1.00 


1.00 


11.523843 


11.455969 


0.19 


0.67 


7T 

on 
20 


1.00 


1.00 


1.00 


11.523843 


11.386608 


0.19 


0.69 


7T 

on 
20 


1.00 


1.00 


1.00 


11.523843 


11.344747 


0.61 


0.23 


7T 


1.00 


1.00 


1.00 


11.523843 


11.516685 


0.63 


0.19 


7T 


1.00 


1.00 


1.00 


11.523843 


11.489791 


0.63 


0.21 


7T 


1.00 


1.00 


1.00 


11.523843 


11.436105 


0.65 


0.17 


7T 


1.00 


1.00 


1.00 


11.523843 


11.461786 


0.65 


0.19 


7T 


1.00 


1.00 


1.00 


11.523843 


11.392115 


0.65 


0.21 


7T 


1.00 


1.00 


1.00 


11.523843 


11.349366 


0.67 


0.15 


7T 


1.00 


1.00 


1.00 


11.523843 


11.472422 


0.67 


0.17 


7T 


1.00 


1.00 


1.00 


11.523843 


11.383978 


0.67 


0.19 


7T 


1.00 


1.00 


1.00 


11.523843 


11.324970 


0.67 


0.21 


7T 


1.00 


1.00 


1.00 


11.523843 


11.291915 


0.69 


0.13 


7T 


1.00 


1.00 


1.00 


11.523843 


11.522980 


0.69 


0.15 


7T 


1.00 


1.00 


1.00 


11.523843 


11.412094 


0.69 


0.17 


7T 


1.00 


1.00 


1.00 


11.523843 


11.334189 


0.69 


0.19 


7T 


1.00 


1.00 


1.00 


11.523843 


11.284714 



For the parameters a,b,6 given in Table 1 and Table 2, we all have ofo < d\. This completes 
the Proof of Theorem 1.1. □ 

Proof of Theorem 1.2 To get Theorem 1.2, we are going to find a test loop q G A + such 
that f(q) < c?3. Then the minimizer of / on A + must be a noncollision solution if d% < d\. 



Let a > 0, 6 G [0, 2vr) and 

q — q\ = ae x 



-i(^t+fl) 



(3.25) 



Hence 



Q - Q2 



qi + ae ^i(^t+e) 



rie v 



<?2 



+ ae y 



(3.26) 



Q-Q3 = qi + ae^T^ 6 



rie v 



93 



(3.27) 



It is easy to see that q G A + and 



\q - ?i I 



/2vr\2 



qi\ 



(3.28) 
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\Q ~ 92 1 



2vr\2 



T J 



) [a 2 + r? + rf - 2rir 2 co,s(6l 2 - 6>i) + 2aricos(6»i - 0) - 2ar 2 cos{9 2 - 9)] , (3.29) 



k - 92 1 = [a 2 + f\ + Hj — 2rir 2 cos(0 2 - #i) + 2aricos( 



2ar 2 cos(6» 2 - 6»)] 2 , (3.30) 



2ir\ 2 

[ —J [a 2 + r 2 + r 2 - 2r 1 r 3 cos(0 3 - X ) + 2ancos(0i - 0) - 2ar 3 cos(^3 - 0)] , (3.31) 



|g - ?3 | = [a 2 + r\ + r| - 2rir 3 cos(6> 3 - 0i) + 2aricos(0i - 9) - 2ar 3 cos(6 3 - 9)] 

, 2 



k'il 2 = (y) r i' l^l 2 = (y) r2 ' '^ 3 ' 2 = V 



= /2vry 



'3- 



cos (63 — 6 1 



Therefore by ([338]) - (13331 . we get 

1 /" T ^ r 1 , ■ ■ ,2 m 1, 1 /" T A . . .0 , 

/(9) = "TT / > TO; — O — (7j + 1 r TT / > ™? 9i dt 

M Jo {ri ^ |g-?i|J 2m; ^ 

27r 2 r 2 m 2 + m 3 - mi 2 2m 2 r 2 cas(0 2 — #i) + 2m 3 r 3 

2(m 2 + m 3 ) .„ n . 2m,2r2Cos{92 - 9) + 2m 3 r 3 cos 
+ ancosiO, -9) 

H h m 2 / [a 2 + r 2 + rf — 2rir2Cos(9 2 — $i) + 2ar\cos(9i — 

a Jo 

_i /-T 

-2ar 2 cos(0 2 - 0)] 2 dt + m 3 / [a 2 + r 2 + r| - 2rir 3 cos(0 3 - 9i) 

Jo 

+2ar l cos{9 l -9)- 2ar 3 cos(9 3 - 0)]~*dt 
= d 3 (a,9). 



(3.32) 
(3.33) 



■n 



(03 ~ g) ■ 
a 



(3.34) 



In order to estimate d 3 , we have computed the numerical values of d 3 = f(q) over some selected 
test loops. The computation of the integral that appears in (|3.34p has been done using the function 
{quad} of Mathematica 7.1 with an error less than 10 -6 . Let T = 1, the results of the numerical 
explorations are given in Table 3 with M = 1 and Table 4 with mi = m 2 = m 3 = 1. 
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Table 3: Parameters for test loops for Theoerm 1.2 



a 


e 


ml 


m2 


m3 


dl 


d3 


0.17 


n 
z 


0.10 


0.75 


0.15 


5.062791 


5.060773 


0.17 


7T 
o 

z 


0.10 


0.77 


0.13 


5.083903 


5.071551 


0.17 


7T 
o 

z 


0.10 


0.78 


0.12 


5.094969 


5.077450 


0.17 


7T 
o 

z 


0.10 


0.80 


0.10 


5.118123 


5.090270 


0.17 


7T 

2 


0.15 


0.53 


0.32 


5.051742 


5.050040 


0.17 


7T 

2 


0.15 


0.57 


0.28 


5.068768 


5.046398 


0.17 


7T 
o 

z 


0.15 


0.60 


0.25 


5.085112 


5.047242 


0.17 


7T 
o 

z 


0.15 


0.65 


0.20 


5.119162 


5.055458 


0.17 


7T 
o 

z 


0.15 


0.70 


0.15 


5.161725 


5.072186 


0.17 


7T 
o 

z 


0.15 


0.72 


0.13 


5.121130 


5.081261 


0.17 


7T 
o 

z 


0.20 


0.31 


0.49 


5.176554 


5.175168 


0.17 


7T 
o 

z 


0.20 


0.35 


0.45 


5.167020 


5.144967 


0.17 


7T 
o 

z 


0.20 


0.40 


0.40 


5.162763 


5.114876 


0.17 


7T 

2 


0.20 


0.50 


0.30 


5.179789 


5.080232 


0.17 


7T 
o 

z 


0.20 


0.55 


0.25 


5.201070 


5.075680 


0.17 


7T 
o 

z 


0.20 


0.60 


0.20 


5.230864 


5.079639 


0.19 


7T 
o 

z 


0.25 


0.22 


0.53 


5.249837 


5.237465 


0.19 


7T 
o 

z 


0.25 


0.25 


0.50 


5.325541 


5.202291 


0.19 


7T 
o 

z 


0.25 


0.30 


0.45 


5.308516 


5.150479 


0.19 


7T 
o 

z 


0.25 


0.35 


0.40 


5.300003 


5.107178 


0.19 


7T 
o 

z 


0.25 


0.62 


0.13 


5.041112 


5.020454 


0.19 


7T 
o 

z 


0.30 


0.22 


0.48 


5.230258 


5.222385 


0.19 


7T 

2 


0.30 


0.25 


0.45 


5.308516 


5.189765 


0.19 


7T 
o 

z 


0.30 


0.30 


0.40 


5.445755 


5.142208 


0.19 


7T 
o 

z 


0.30 


0.35 


0.35 


5.441499 


5.103164 


0.19 


7T 
o 

z 


0.30 


0.56 


0.14 


5.036553 


5.032137 


0.21 


7T 
o 

z 


0.35 


0.21 


0.44 


5.192935 


5.184596 


0.21 


7T 
o 

z 


0.35 


0.29 


0.36 


5.412519 


5.092092 


0.21 


7T 
o 

z 


0.35 


0.39 


0.26 


5.327621 


5.007107 


0.21 


7T 
o 

z 


0.35 


0.48 


0.17 


5.091316 


4.959734 


0.21 


7T 

2 


0.35 


0.53 


0.12 


4.971952 


4.945333 


0.21 


7T 

3 


0.40 


0.28 


0.32 


5.386433 


5.342981 


0.21 


7T 

Q 
■1 


0.40 


0.32 


0.28 


5.386433 


5.287294 


0.21 


7T 

Q 


0.40 


0.36 


0.24 


5.271874 


5.237055 


0.21 


7T 

Q 
•1 


0.40 


0.38 


0.22 


5.216638 


5.213978 


0.23 


7T 

o 
Z 


0.45 


0.19 


0.36 


5.139742 


5.127834 


0.23 


7T 

2 


0.45 


0.29 


0.26 


5.337836 


5.003006 


0.23 


7T 

2 


0.45 


0.37 


0.18 


5.112805 


4.927660 


0.23 


7T 

2 


0.45 


0.46 


0.09 


4.885693 


4.868944 


0.23 


7T 

2 


0.50 


0.18 


0.32 


5.123871 


5.108878 


0.23 


7T 

2 


0.50 


0.23 


0.27 


5.266218 


5.044762 


0.23 


7T 

2 


0.50 


0.29 


0.21 


5.208258 


4.979058 


0.23 


7T 

2 


0.50 


0.37 


0.13 


4.990036 


4.910522 


0.23 


7T 

2 


0.50 


0.41 


0.09 


4.889098 


4.884426 
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Table 4: Parameters for test loops for Theoerm 1.2 



a 


e 


ml 


m2 


m3 


dl 


d3 


0.21 


7T 

o 
z 


1.00 


1.00 


1.00 


11.523843 


11.327950 


0.23 


7T 

o 
z 


1.00 


1.00 


1.00 


11.523843 


11.036769 


0.23 


7T 

Q 
"1 


1.00 


1.00 


1.00 


11.523843 


11.336568 


0.25 


7T 

o 
z 


1.00 


1.00 


1.00 


11.523843 


10.821272 


0.25 


7T 


1.00 


1.00 


1.00 


11.523843 


11.187475 


0.25 


7T 
4 


1.00 


1.00 


1.00 


11.523843 


11.453195 


0.27 


7T 

o 
z 


1.00 


1.00 


1.00 


11.523843 


10.667031 


0.27 


7T 


1.00 


1.00 


1.00 


11.523843 


11.107374 


0.27 


7T 
4 


1.00 


1.00 


1.00 


11.523843 


11.411685 


0.29 


7T 

o 
z 


1.00 


1.00 


1.00 


11.523843 


10.563849 


0.29 


7T 

Q 


1.00 


1.00 


1.00 


11.523843 


11.085761 


0.29 


7T 
A 

4 


1.00 


1.00 


1.00 


11.523843 


11.430090 


0.31 


7T 

o 
z 


1.00 


1.00 


1.00 


11.523843 


10.504424 


0.31 


7T 

Q 
■J 


1.00 


1.00 


1.00 


11.523843 


11.114860 


0.31 


7T 
4 


1.00 


1.00 


1.00 


11.523843 


11.500414 


0.33 


7T 

2 


1.00 


1.00 


1.00 


11.523843 


10.483477 


0.33 


7T 

3 


1.00 


1.00 


1.00 


11.523843 


11.188786 


0.35 


7T 

2 


1.00 


1.00 


1.00 


11.523843 


10.497161 


0.35 


7T 

3 


1.00 


1.00 


1.00 


11.523843 


11.302997 


0.37 


7T 

2 


1.00 


1.00 


1.00 


11.523843 


10.542652 


0.37 


7T 

3 


1.00 


1.00 


1.00 


11.523843 


11.453926 


0.39 


7T 

2 


1.00 


1.00 


1.00 


11.523843 


10.617860 



For the parameters a, 9 given in Table 3 and Table 4, we all have 0(3 < d\. This completes the 
Proof of Theorem 1.2. □ 
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